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Abstract Topology in momentum space is the main characteristics of the ground 
states of a system at zero temperature, the quantum vacua. The gaplessness of 
fermions in bulk, on the surface or inside the vortex core is protected by topology, 
' and is not sensitive to details of the microscopic physics (atomic or trans-Planckian). 

Irrespective of the deformation of the parameters of the microscopic theory, the en- 
^ I ergy spectrum of these fermions remains strictly gapless. This solves the main hier- 

. archy problem in particle physics: for fermionic vacua with Fermi points the masses 

of elementary particles are naturally small. The quantum vacuum of Standard Model 
is one of the representatives of topological matter alongside with topological super- 
fluids and superconductors, topological insulators and semi-metals, etc. There is a 
number of of topological invariants in momentum space of different dimensions. 
They determine universality classes of the topological matter and the type of the 
effective theory which emerges at low energy. In many cases they also give rise to 
emergent symmetries, including the effective Lorentz invariance, and emergent phe- 
nomena such as effective gauge and gravitational fields. The topological invariants 
in extended momentum and coordinate space determine the bulk-surface and bulk- 
vortex correspondence. They connect the momentum space topology in bulk with 
the real space. These invariants determine the gapless fermions living on the surface 
of a system or in the core of topological defects (vortices, strings, domain walls, 
solitons, monopoles, etc.). The momentum space topology gives some lessons for 
quantum gravity. In effective gravity emerging at low energy, the collective variables 
are the tetrad field and spin connections, while the metric is the composite object of 
tetrad field. This suggests that the Einstein-Cartan-Sciama-Kibble theory with tor- 
sion field is more relevant. There are also several scenarios of Lorentz invariance 
violation governed by topology, including splitting of Fermi point and development 
of the Dirac points with quadratic and cubic spectrum. The latter leads to the natural 
emergence of the Hof ava-Lifshitz gravity. 
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1 Introduction 

Topological approach seeks to recover the Standard Model of particle physics and 
general relativity as the low-energy phenomena emerging in the complicated micro- 
scopic media called the quantum vacuum. Momentum space topology distributes 
the quantum vacua into universality classes, which are characterized by different 
types of topological invariants, expressed in terms of propagators - Green's func- 
tions. The similar invariants characterize the ground states of topological materials: 
topological superfluids, superconductors, insulators, semi-metals, etc. The only dif- 
ference between the topological materials and the quantum vacuum of the Standard 
Model is that the latter obeys several symmetries, which are absent in condensed 
matter physics. These are the Lorentz invariance and different types of gauge sym- 
metry. However, in principle, the symmetries of the Standard Model can arise as 
the emergent phenomena. This is precisely what happens in the universality class of 
vacua, which are characterized by the existence of the topologically protected Fermi 
point in momentum space - the point where the energy of the fermionic excitations 
is necessarily nullified due to non-zero topological invariant. Such point exists in 
the vacuum of the Standard Model in the symmetric phase above the electroweak 
interaction, and also in some topological materials, whose ground states also be- 
long to the Fermi point universality class. These are superfluid ^He in the phase A, 
and topological semimetals. That is why we call the symmetric phase of the Stan- 
dard Model as quantum vacuum in the semi-metal state. At low energy, i.e. in the 
vicinity of the point node in the spectrum of fermionic excitations, these excitations 
behave as relativistic particles - chiral (left handed or right handed) Weyl fermions 
- interacting with effective gauge fields and gravitational fields. In principle, all the 
ingredients of Standard Model and general relativity can be recovered as phenomena 
naturally emerging in the background of the vacuum (ground state), which belongs 
to the Fermi point universaUty class. 

Here we consider also the other universaUty classes of topological media, illus- 
trating them on condensed matter and particle physics examples. Among them is the 
vacuum of Standard Model below the electroweak transition, where the momentum- 
space topology prescribes the appearance of masses for all fermionic excitations: 
quarks and leptons. This massive phase of the Standard Model vacuum is also topol- 
gically nontrivial: it belongs to universality class, which includes such topological 
media as superfluid ^He in the phase B and three-dimensional topological insula- 
tors obeying time reversal symmetry. The fermionic spectrum in all these systems is 
fully gapped, but there are topologically protected gapless (massless) states on the 
surface of these media or at the interfaces. That is why we call the massive phase of 
the Standard Model as the quantum vacuum in the insulating state. 

Another interesting class of topological media includes the vacua where the 
fermionic excitations are exotic gapless Dirac particles with non-relativistic quadratic, 
cubic or higher order spectrum. The effective quantum electrodynamics, which 
emerges in such vacua, experiences the so-called anisotropic scaUng, in which the 
space and time transform in different ways. The topological mechanism of emergent 
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anisotropic scaling can serve as a source of the recently suggested quantum gravity 
at short distances (the so-called Hofava-Lifshitz gravity). 

Some sections are devoted to the fermionic bound states leaving on the bound- 
ary of the topological media, on the interfaces between the vacua with different 
topological invariants, or inside the topological defects, such as strings in relativis- 
tic theories and vortices in superfluids and superconductors. The topology can say 
may or may not these bound states have exactly zero energy. In superfluids and 
supercondictors such zero-energy states behave as Majorana fermions, which are 
still elusive in particle physics. The topological consideration of bound states sug- 
gests another scenario of emergence of the chiral fermions of the Standard Model: 
they may emerge as the fermion zero modes at the interface between two different 
insulating quantum vacua in the 4-+1 space-time. 



2 Quantum vacuum as topological medium 

2.1 Symmetry vs topology 

There is a fundamental interplay of symmetry and topology in physics, both in con- 
densed matter and relativistic quantum fields. Traditionally the main role was played 
by symmetry: gauge symmetry of Standard Model and GUT; symmetry classifica- 
tion of condensed matter systems such as solid and liquid crystals, magnets, super- 
conductors and superfluids; universality classes of spontaneously broken symmetry 
phase transitions; etc. The last decades demonstrated the opposite tendency in which 
topology is becoming primary being the main characteristics of quantum vacua in 
relativistic quantum fields and their condensed-matter counterparts - ground states 
of the condensed-matter systems at T = 0, see reviews |[Tl|2l[3]|4l and earlier papers 

13 |6l El 111 IE [13 III] El [!3 [II [B] [111 [nl [11] . 

Topology describes the properties of a system, which are insensitive to the details 
of the microscopic physics (atomic physics in condensed matter and Planck-scale 
physics of quantum vacuum). The momentum-space topology, which we discuss in 
this Chapter, determines universality classes of the topological media and the ef- 
fective quantum field theory which emerges in such media at low energy and low 
temperature, including the type of the energy spectrum of fermionic excitations (in 
condensed matter systems fermionic excitations play the role of elementary parti- 
cles). The topology also gives rise to emergent symmetry, the symmetry which did 
not exist on a microscopic level. Examples are provided by the point nodes (zeroes) 
in the energy spectrum of fermionic excitations: some of the nodes are protected 
by topology, i.e. they are robust to the deformations of the system. Close to such 
nodes the spectrum of fermions becomes "relativistic", i.e. the spectrum forms the 
Dirac cone, and fermionic excitations behave as Weyl, Dirac or Majorana particles. 
As a result the effective relativistic quantum field theory emerges, in which bosonic 
collective modes give rise to effective gauge field and effective metric. All this is 
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the consequence of the topological theorem - the Atiyah-Bott-Shapiro construction 

im. 

Among the existing and potential representatives of topological materials one can 
find those in which the spectrum of fermionic excitations has a gap (the gap corre- 
sponds to mass of a relativistic particle). These systems include 3D topological band 
insulators |3l; fully gapped superluid ^He-B lfT4l [T9l l20l : 2D materials exhibiting 
intrinsic (i.e. without external magnetic field) quantum Hall and spin-Hall effects, 
such as gapped graphene ifTsl : thin film of superluid ^He-A and quasi 2D planar 
phase of triplet superfluid lfT6l 1211 [iTl ; and chiral superconductor Sr2Ru04 ll22ll . 
These materials have the topological properties similar to that of the quantum vac- 
uum of Standard Model in its massive phase, i.e. in the state below the electroweak 
transition, where all elementary particles are massive |[23l. The latter state can be 
called the "insulating" state of the relativistic quantum vacuum. 

The gapless topological media are represented by superluid ^He-A HI; topo- 
logical semimetals EH |25] EH |27| ; gapless graphene 118] [lH 12 EO] EH ; nodal 
cuprate superconductors ll28l : and noncentrosymmetric |[32l, l33l superconductors. 
These materials are similar to the quantum vacuum of Standard Model in the state 
above the electroweak transition, where all elementary particles are massless lfn i23l . 
This state can be called the "semimetal" phase of the relativistic quantum vacuum. 



2.2 Green's function vs order parameter 

Topology operates in particular with integer numbers - topological invariants or 
topological charges - which do not change under small deformation of the system. 
The topological invariants appear in very different situations. The commonly known 
are topological charges, which describe the topological defects - the inhomogeneous 
objects such as quantized vortices, dislocations, domain walls, hedgehogs, solitons, 
etc. in condensed matter systems or the cosmic strings, magnetic monopoles, instan- 
tons, etc. in particle physics. The topological charge protects the topological defects 
from destruction, while the summation law of the topological charges regulates pro- 
cesses of scattering, merging, splitting, and other transformations of these defects in 
their dynamics. 

Here we consider the topological invariants of different type: they describe the 
homogeneous ground states or the homogeneous quantum vacua. The conservation 
of these topological charges in particular protects from destruction different types of 
nodes in the fermionic energy spectrum of a system. The node in the spectrum is an 
object in momentum space, and thus we call these invariants as momentum-space 
invariants, as distinct from the real-space invariants describing topological defects. 
In other words, the real-space invariants describes the topologically nontrivial inho- 
mogeneous configurations of the order parameter fields f (r,f) in space-time, while 
the momentum-space invariants describe the nontrivial momentum-space config- 
uration of the propagator - Green's function G(p, o) or other response function, 
which characterizes the homogeneous ground state of a system (the vacuum state) 
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Khodel-Shaginyan flat band: Tt-vortex in p-space 

Fig. 1 Topological matter, represented in terins of topological objects in moinentum space, (top 
lefty. Fermi surface is the momentum-space analog of the vortex line: the phase of the Green's 
function changes by Itc around the element of the line in (co,p)-space. (top middle): Fermi point 
(Weyl point) is the counterpart of a hedgehog and a magnetic monopole. The hedgehog in this 
figure has integer topological charge N = +1, and close to this Fermi point the fermionic quasi- 
particles behave as Weyl fermions. Nontrivial topological charges in terms of Green's functions 
support the stability of the Fenni surfaces and Weyl points with respect to perturbations including 
interactions (11 1341 . In terms of the Berry phase 1361 the Fermi point represents the p-space coun- 
teipart of Dirac magnetic monopole with unobservable Dirac string (see Ref . 1371 and Fig. 1 1 .4 in 
fTl). (top right): Topological insulators and fully gapped topological superfluids/superconductors 
are textures in momentum space: they have no singularities in the Green's function and thus no 
nodes in the energy spectrum in the bulk. This figure shows a skyrmion in the two-dimensional 
momentum space, which characterizes two-dimensional topological insulators exhibiting intrinsic 
quantum Hall or spin-Hall effect, (bottom left): Flat band emerging in strongly interacting systems 
1411 . This dispersionless Fermi band is analogous to a soliton terminated by half-quantum vortices: 
the phase of the Green's function changes by n around the edge of the flat band 1421 . (bottom right): 
Fermi arc on the surface of ^^He-A 1431 and of topological semi-metals with Weyl points |27||26l 
and flat band inside the vortex core of^He-Al44l serve as the momentum-space analog of a Dirac 
string terminating on a monopole. The Fenni surface formed by the surface bound states terminates 
on the points where the spectrum of zero energy states merge with the continuous spectrum in the 
bulk, i.e. with the Weyl points. 
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im [T8l |34] [35l . The Green's function is generally a matrix with spin indices. In ad- 
dition, it may have the band indices (in the case of electrons in the periodic potential 
of crystals). 

Due to topological stability, the nodes in the spectrum survive when the interac- 
tion between the fermions is introduced and/or modified. Thus the momentum-space 
topology is the main reason why there are gapless quasiparticles in condensed mat- 
ter and (nearly) massless elementary particles in our Universe. Topology explains 
the masslessness of the quantum vacuum in its semimetal phase, and allows to solve 
the hierarchy problem in particle physics: if the masslessness (or gaplessness) of the 
quantum vacuum is not protected, the natural values of the masses of elementary 
particles are on the order of the Planck energy scale, which is in huge disagreement 
with masses of elementary particles in the insulating vacuum. However, if masses 
appear as a result of the spontaneously broken symmetry in the initially gapless 
vacuum, they can be small. 

The momentum-space topological invariants are in many respects similar to their 
real-space counterparts, which describe topological defects in condensed matter sys- 
tems and quantum vacuum (see Fig.[T]i. In particular, the point node in the spectrum 
(the Fermi point or Weyl point, see Sec. [3]i is the counterpart of the real-space point 
defects, such as hedgehog in ferromagnets or magnetic monopole in particle physics 
(Fig.[T]fo/7 center). The summation law of the momentum-space topological charges 
regulates processes of scattering, merging, splitting, and other transformations of the 
nodes in spectrum during deformation of the vacuum state. Another example is the 
Fermi surface in metals, which we consider in Sec. 12.31 It is topologically stable, 
because it is analogous to the vortex loop in superfluids or superconductors (Fig.[T] 
top left). 

The fully gapped topological matter, such as topological insulators and fully 
gapped topological superfluids have no nodes in the spectrum and thus no singu- 
larities in the Green's function. They are the counterpart of the real-space objects 
which do not have singularities, but still are topologically non-trivial - the so-called 
textures or skyrmions. Thus the topological insulator and the vacuum of Standard 
Model in insulating phase can be considered as the momentum-space skyrmion 
(Fig. 1 top right). 

The further extension of topology to the combined real + momentum space al- 
lows us to consider topologically protected spectrum of the so called fermion zero 
modes - fermions living on the real-space topological objects such as domain walls, 
strings and monopoles. The fermion zero modes are protected by topological in- 
variants expressed in terms of the Green' s function in the extended phase space, 
G(co,p;f,r) ||38llIl|39l[I]|40l|3l. 



2.3 Fermi surface as topological object 

We shall start with gapless vacua - the vacua whose fermionic excitations has no 
gap. In the gapless media the fermionic degrees of freedom are not frozen out at 
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Fig. 2 Nodes of co-dimension 1 in 2+1 systems. Green's function has singularities on line £0 = 0, 
P? + p^, = in the three-dimensional space ((B,p, , p, ). Stability of Fermi surface is protected by 
the invariant {TJ which is represented by integral over an arbitrary contour C around the Green's 
function singularity. This is applicable to nodes of co-dimension 1 in any D+\ dimension. For 
Z) = 3 the nodes form conventional Fenni surface in metals and in normal ^^He. 

low temperature. This is very important for practical purposes in electronic devices 
and for the existence of life in our Universe. It is the topology which provides the 
protection for the gaplessness in condensed matter systems and masslessness of 
elementary particles in Standard Model. 

For the topological classification of the gapless vacua, the Green's function is 
considered on imaginary frequency axis. On the real axis singularities are always 
present, for example in the simplest case of non-interacting fermions the poles in 
the Green's function, G = 1/ (o — e(p)), are present even in the fully gapped sys- 
tem. To distinguish between the gapped and gapless states of the vacuum, and also 
between different types of gapless vacua, we should move to the imaginary axis. 
Then the simplest Green's function becomes G = 1/ {i(0 — e(p)), which may have 
singularity at O) = and e(p) = 0. The latter is possible only for the gapless vacua, 
where the energy of excitations e(p) has zeroes. Thus zeroes in the spectrum of 
fermionic excitations are manifested as singularities in the Green's function, and we 
are interested in such singularities in the Green's function, which are protected by 
topology, i.e. are described by topological charge and thus are robust to perturba- 
tions of the vacuum state. The transition to imaginary axis thus allows us to consider 
only the relevant singularities in the Green's function and to avoid the trivial singu- 
larities on the mass shell. 

Zeroes in the energy spectrum of fermionic excitations may form different p- 
dimensional manifolds in momentum space: there can be zero-dimensional points 
p — Q, one-dimensional lines of zeroes p = I, two-dimensional surfaces p ~2 and 
even the whole three-dimensional bands with zero energy p ^3. Since the dimen- 
sion D of the momentum space can be also different, for topological classification it 



8 



G.E. Volovik 



is more instructive to use the co-dimension of zeroes. By co-dimension we denote 
D — p, i.e. the dimension D of p-space minus dimension p of the nodes in energy 
spectrum. 

We start with zeroes of co-dimension 1 . This refers to two-dimensional Fermi 
surface in three-dimensional metal (3 — 2 = 1), D=l Fermi line in D=2 systems 
(2—1 = 1) and Fermi point in D = 1 systems (1 — = 1). The general analysis ifTSl 
demonstrates that topologically stable nodes of co-dimension 1 are described by the 
group Z of integers. The corresponding winding number is expressed analytically 
in terms of the Green's function |[T]: 

N = tr <f^^.Gico,p)d,G-\co,p) . (1) 

Here the integral is taken over an arbitrary contour C around the Green's function 
singularity in the D+l momentum-frequency space. See Fig.|2]for D = 2. Example 
of the Green's function in any dimension D is scalar function G^'(a),p) = i(0 — 
vf(|p| — Pf)- For D = 2, the singularity with winding number = 1 is on the line 
(0 = 0, p^ + p^, = p^p, which represents the one-dimensional Fermi surface. 

Due to nontrivial topological invariant, Fermi surface survives the perturbative 
interaction and exists even in systems without poles in the Green's function and 
quasiparticles are not well defined. The systems without poles include the so-called 
marginal and Luttinger liquids in condensed matter and the so-called unparticles in 
relativistic theories (see Sec. |3.2| l. 



3 Vacuum in a semi-metal state 

Here we consider the topology of the quantum vacuum of Standard Model in its 
symmetric phase above the electroweak transition, where all fermions are massless 
Weyl fermions. Standard Model vacuum obeys the relativistic invariance and thus 
cannot belong to the universality class of vacua with Fermi surface: the latter vio- 
lates the Lorentz symmetry. The gapless (massless) phase of Standard Model can 
be consistent with Lorentz symmetry only if the vacuum of Standard Model be- 
longs to the universality class of vacua, which is characterized by point nodes in 
the fermionic spectrum - the Fermi points. The Fermi point has the co-dimension 
3. The Fermi point can be of the chiral type, when the fermionic excitations near 
the Fermi point behave as relativistic left-handed or right-handed Weyl fermions, 
such as left neutrino, left electron, left quarks and their right-handed counterparts 
P31 [l1. In this case the Fermi point is called the Weyl point. The Fermi point can 
be also of the Dirac type, when the fermionic excitations behave as massless Dirac 
fermions. There is also the class of quantum vacua in which the nodal point obeys 
the Z2 topology, with the topological charge obeying the summation rule 1 + 1=0. 
In such vacuum the fermionic excitations behave as relativistic massless Majorana 
particles ifTsl . 
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Fermi point is the Green's function singularity described by the following topolog- 
ical invariant expressed via integer valued integral over the surface a around the 
singular point in the 4-momentum space = (fi),p) HI: 

^ = ^ I ^^Pp G ' ^^P^ G ' ^^/^v G '. (2) 

If the invaraint (O is nonzero, the Green's function has a singularity inside the sur- 
face C7, and this means that fermions are gapless. The typical singularities have 
topological charge A^=:+l orA^=— l.In the vicinity of the Fermi point with such 
topological invariant the spectrum of fermionic excitations reproduces the spectrum 
the chiral Weyl fermions, right-handed or left-handed respectively. This is the con- 
sequence of the so-called Atiyah-Bott-Shapiro construction ifTSl . which leads to the 
following general form of expansion of the inverse fermionic propagator near the 
Fermi point with A^ = +lorA'=— 1: 

G-\p,)=e^ar"{pp-pf^) + ---. (3) 

Here = (1 , CJj, CJy, o;) are Pauli matrices (or Dirac matrices in the more general 
case); the expansion parameters are the vector p^^ ^ indicating the position of the 
Fermi point in momentum space where the Green's function has a singularity, and 
the matrix ; ellipsis denote higher order terms in expansion. The expansion results 
in the linear spectrum of fermionic excitations, which induces the effective Lorentz 
invariance. 



3.2 Emergent relativistic fermionic matter 

By continuously adiabatic deformation one may transform equation (O to the equa- 
tion which describes the relativistic Weyl fermions 

G-\p^) = i(0+Na-p+--- , N = ±l. (4) 

Here the position of the Fermi point is shifted to p^^^ = 0; the matrix is deformed 
to unit matrix; and ellipsis denote higher order terms in (O and p. This means that 
close to the Fermi point with = + 1 , the low energy fermionic excitations behave 
as right handed relativistic particles, while the Fermi point with = — 1 gives rise 
to the left handed particles. They obey the following effective Weyl Hamiltonian 



(5) 
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Fig. 3 Weyl point as the conical point of level crossing. The negative energy branch occupied by 
the fermionic excitations touches the positive energy branch. 



The spectrum of this Hamihonian is shown in Fig. [3], which demonstrates that 
the Fermi or Weyl point is an example of the exceptional point of level crossing 
analyzed by von Neumann and Wigner P6l . This analysis demonstrates that two 
branches of spectrum, which have the same symmetry, may touch each other at the 
conical (or diabolical) point in the three-dimensional space of parameters. In other 
words the twofold degeneracy has co-dimension = 3. In our case the two branches 
correspond to positive and negative energy states, while the three parameters are the 
components of momentum p,, py and |[8]|32l. Topological invariants for points 
at which the branches of spectrum merge were introduced by Novikov 1471 . 

However, this is not the whole story. Topological invariant of the Weyl point does 
not change if the Green's function is multiplied by any power of p^ + o^, which 
violates the simple pole structure of the propagator in Eq.(|4li. Such violation may 
occur due to the infrared divergences, which take place in quantum field theories. In 
this case in the vicinity of Fermi point one has 



with 7 7^ 1 . In particle physics, the hypothetical particles described by the modified 
Green's function with 77^ 1 are called unparticles 1481 |49l . This modification does 
not change the topology of the propagator: the topological charge of singularity is 
for arbitrary parameter 7 lO . For fermionic unparticles one has 7 = 5 /2 — t/y, where 
du is the scale dimension of the quantum field. 

The main property of the vacua with Weyl points is that according to (|4|i, close 
to the Weyl points the massless relativistic fermions emerge. This is consistent with 
the fermionic content of our Universe, where all the elementary particles - left- 
handed and right-handed quarks and leptons - are Weyl fermions. Such a coinci- 
dence demonstrates that the vacuum of Standard Model is the topological medium 
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of special universality - the class of vacua which have topologically protected Fermi 
points. This solves the hierarchy problem, since the value of the masses of elemen- 
tary particles in the vacua of this universality class is identically zero. 

Let us suppose for a moment, that there is no topological invariant which protects 
massless fermions, and the vacuum of our Universe belongs to the class of the fully 
gappedvacua. Then the natural masses of fermions must be on the order of Planck 
energy scale: M ^Ep ^ lO'^ GeV. In such a natural Universe, where all masses are 
of order E-p, all fermionic degrees of freedom are completely frozen out because of 
the Bolzmann factor e^'^l'^ , which is about at the temperature corresponding 

to the highest energy reached in accelerators. There is no fermionic matter in such 
a Universe at low energy. That we survive in our Universe is not the result of the 
anthropic principle (the latter chooses the Universes which are fine-tuned for life 
but have an extremely low probability). Our Universe is also natural and its vacuum 
is generic, but it belongs to the universality class of vacua with Fermi points. In 
such vacua the masslessness of fermions is protected by topology (combined with 
symmetry, see below). 



3.3 Emergent gauge fields 

The vacua with Fermi-point suggest a particular mechanism for emergent symme- 
try. The Lorentz symmetry is simply the result of the linear expansion: this sym- 
metry becomes better and better when the Fermi point is approached and the non- 
relativistic higher order terms in Eq.© may be neglected. This expansion demon- 
strates the emergence of the relativistic spin, which is described by the Pauli matri- 
ces. It also demonstrates how gauge fields and gravity emerge together with chiral 
fermions. The expansion parameters z?^"' and may depend on the space and time 
coordinates and they actually represent collective dynamic bosonic fields in the vac- 
uum with Fermi point. The vector field p^*'^ in the expansion plays the role of the 
effective U{\) gauge field acting on fermions. 

For the more complicated Fermi points with |A^| > 1 the shift p^'^ becomes the 
matrix field; it gives rise to effective non-Abelian (Yang-Mills) SU {N) gauge fields 
emerging in the vicinity of Fermi point, i.e. at low energy IH. For example, the 
Fermi point with N = 2 may give rise to the effective SU{2) gauge field in addition 
to the effective U{1) gauge field 

G'\Pn)=ear" {pji -giAji -g2A|3 -t) + higher order terms , (7) 
where t are Pauli matrices corresponding to the emergent isotopic spin. 
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3.4 Emergent gravity 

The matrix field in (|7|i acts on the (quasi)particles as the field of vierbein, and thus 
describes the emergent dynamical gravity field. As a result, close to the Fermi point, 
matter fields (all ingredients of Standard Model: chiral fermions and Abelian and 
non-Abelian gauge fields) emerge together with geometry, relativistic spin, Dirac 
matrices, and physical laws: Lorentz and gauge invariance, equivalence principle, 
etc. This is the result of the natural coarse graining: in the low energy corner, when 
the Fermi point is approached, the huge number of microscopic (Planckian) degrees 
of freedom of the physical quantum vacuum is reduced to rather few low-energy 
degrees of freedom (chiral fermions, gauge fields and gravity field). All other de- 
grees of freedom are highly massive and are completely frozen out at low energy 
and temperature. 

In such vacua, gravity emerges together with matter If this Fermi point mech- 
anism of emergence of physical laws works for our Universe, then the so-called 
"quantum gravity" does not exist. The gravitational degrees of freedom become 
separated from all other degrees of freedom of quantum vacuum only at low energy. 

In this scenario, classical gravity is a natural macroscopic phenomenon emerging 
in the low-energy corner of the microscopic quantum vacuum, i.e. it is a typical and 
actually inevitable consequence of the coarse graining procedure discussed above. 
It is possible to quantize gravitational waves to obtain their quanta - gravitons, since 
in the low energy corner the results of microscopic and effective theories coincide. 
It is also possible to obtain some (but not all) quantum corrections to Einstein equa- 
tion and to extend classical gravity to the semiclassical level. But one cannot obtain 
"quantum gravity" (i.e. microscopic physics of the vacuum) by quantization of Ein- 
stein equations for the gravitational degrees of freedom, since all other degrees of 
freedom of the quantum vacuum will be missed in this procedure. 



3.5 Topological invariant protected by symmetry in Standard Model 

Standard Model contains equal number hr ~ ni of right and left Weyl fermions, 8 
right and 8 left in each generation: «r = «l = 8«g, where iig is the number of genera- 
tions (we do not consider Standard Model with Majorana fermions, and assume that 
in the insulating state of Standard Model neutrinos are Dirac fermions). For such 
Standard Model the total topological charge, obtained after trace is taken over all 
fermionic flavors in invariant (|2]i, vanishes, = 8;?;, — 8/1^, = 0. Thus the invariant 
^ cannot protect nodes in the spectrum of Standard Model fermions. This would 
mean that any interaction between the fermions may make them massive. 

However, in the symmetric phase of Standard Modle above the electroweak tran- 
sition this does not happen. This is because there is another topological invariant, 
which takes into account the symmetry of the vacuum. The gapless state of the vac- 
uum with N ^ can be protected by the following topological invariant which is 
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supported by symmetry |[T]: 

(8) 



K I dS" Gdp^ G ' Gdp^ G- ' Gd,,^ G- 



where Kjj is the matrix of symmetry transformation, which either commutes or anti- 
commutes with the Green's function matrix. In the symmetric phase of Standard 
Model there are two relevant symmetries, both are the Z2 groups, ~ \. One 
of them is the center subgroup of SU {2)l gauge group of weak rotations of left 
fermions, where the element K is the gauge rotation by angle 2k, K ~ e'^'^^'-. The 
other one is the group of the hypercharge rotation be angle 6n, K = e'^^^ . In the 
G(224) Pati-Salam extension of the G(213) group of Standard Model, this symme- 
try comes as combination of the Z2 center group of the SU {2)r gauge group for 
right fermions, e"^^3«, and the element e^Tr^B-i) of jjjg ^4 center group of the SU{4) 
color group - the Pm parity (on the importance of the discrete groups in particle 
physics see ifSOl |5T1 and references therein). Each of these two Z2 symmetry op- 
erations changes sign of left spinor, but does not influence the right particles. Thus 
these matrices are diagonal, Kjj = diag(l, 1,..., — 1, — 1,...), with eigen values 1 for 
right fermions and —1 for left fermions. 

In the symmetric phase of Standard Model, both matrices commute with the 
Green's function matrix Gij, as a result Nk in (O is topological invariant: it is robust 
to deformations of Green's function which preserve the symmetry K. The value of 
this invariant A^/f = 16«g, which means that all 16«g fermions are massless. 

Topological invariant protected by symmetry is responsible also for the effect of 
chiral anomaly - the anomalous nonconservation of a chiral current ||52ll53l . This 
anomaly can be the source of the production of baryonic charge B in the early Uni- 
verse, which leads to the excess of matter over anti-matter in the present Universe 
(the so-called electroweak baryogenesis, see the recent review ll54l ). The baryo- 
production can be expressed in terms of the invariant Nk- in case of the creation of 
the baryonic charge by the hypermagnetic U{1) field the relevant symmetry opera- 
tor K in Eq.® is A" = BY^, where Y is the generator of symmetry group. The 
creation rate of baryons in the hypermagnetic U{1) field is 

B=^By-Ey , K^BY\ (9) 

In condensed matter, the chiral anomaly generated by Weyl fermions in a crystal has 
been considered in Ref. Q, while the analog of the baryoproduction in Eq. (|9]l has 
been experimentally tested in the superfluid with Weyl points - the ^He-A | 
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marginal Fermi point 

iV = -1+1 = 




Fermi points witli opposite N marginal Fermi point splits 
annihilate each other & form into two separate 

massive Dirac fermions topologically protected Fermi points 



Fig. 4 (top) In Standard Model the Fermi points with positive N = +i and negative N = — I 
topological charges are at the same point p = 0, fonning the marginal Fermi point with N = 0. 
Symmetry K between the Fermi points prevents their mutual annihilation giving rise to the topo- 
logical invariant (8) with Nk = 2. (bottom left): If symmetry K is violated or spontaneously broken, 
Fermi points annihilate each other and Dirac mass is formed, {bottom right): If Lorentz invariance 
is violated or spontaneously broken, the marginal Fermi point splits |56| . The topological quantum 
phase transition between the state with Dirac mass and the state with splitted Dirac points have 
been observed in cold Fermi gas 1571 . 



3.6 Higgs mechanism vs splitting of Fermi points 

The gapless vacuum of Standard Model is supported by combined action of topol- 
ogy and symmetry K, and also by the Lorentz invariance which keeps all the Fermi 
points at p = 0. 

Explicit violation or spontaneous breaking of one of these symmetries transforms 
the vacuum of the Standard Model into one of the two possible vacua. If, for exam- 
ple, the K symmetry is broken, the invariant (O supported by this symmetry ceases 
to exist, and the Fermi point disappears. All \6n^ fermions become massive (Fig.|4] 
bottom left). This is assumed to happen below the symmetry breaking electroweak 
transition caused by Higgs mechanism where quarks and charged leptons acquire 
the Dirac masses. 

If, on the other hand, the Lorentz symmetry is violated, the marginal Fermi point 
splits into topologically stable Fermi points with non-zero invariant A^, which pro- 
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Fig. 5 Sketch of topologically protected point nodes in momentum space with topological charge 
A' = ±1 in Eq.([2ll in superconfuctors/superfluids of 0{D2) symmetry class |59i . Chiral fermions 
emerge in the vicinity of each Fermi points. They have identical emergent Minkowski metric, 
hut different orientations of dreibein. The simplest realization of dreibein for each of eight chiral 
fermions is shown by aiTows. The vierbein orientations in Ol^Di) symmetry class superconductors 
are more complicated than in Figure: one of the vectors in each vierbein is along the cube main di- 
agonal, see Eq.in). All four pairs of left and right Weyl fermions have the same quantum numbers, 
while their triads can be transformed to each other by rotations and reflection. They are analogous 
to generations of Standard Model fermions. 



tects massless chiral fermions {Vig.^bottom right). Since the invariant does not 
depend on symmetry, the further symmetry breaking cannot destroy the nodes. One 
can speculate that in the Standard Model the latter may happen with the electri- 
cally neutral leptons, the neutrinos ||56ll . Most interestingly, Fermi-point splitting of 
neutrinos may provide a new source of T and CP violation in the leptonic sector, 
which may be relevant for the creation of the observed cosmic matter-antimatter 
asymmetry ll58l . 



3.7 Splitting of Fermi points and problem of generations 

An example of the multiple splitting is provided by the model Hamiltonian for 
fermions in superconfuctors/superfluids in the state which belongs to 0{D2) sym- 
metry class 1591 . The Bogoliubov - de Gennes (BdG) Hamiltonian for fermions in 
this spin-singlet tZ-wave superconductor has the form 

H = i=(/ _ pI)^, + ii2pl-pl-p])x, + -;,2)T2 . (10) 

At < the energy spectrum is fully gapped, for = the node in the spectrum 
appears at p = which at > splits into 8 Fermi points at the vertices of cube in 
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momentum space (see Fig. 13.71 ): 

p(")^Z£(±x±y±z) , „=!,... ,8. (11) 

These nodes have topological charges = ±1 in Eq.(|2]i, and as a result, close to 
each of 8 nodes the Hamiltonian is reduced to the Hamiltonian describing Weyl 
fermions, 4 left and 4 right: 

//W =e<"' • (p-pW)Ti +e("' • (p-pW)T2 + e("^ . (p-pW)T3. (12) 

Each Weyl fermion has its own triad (dreibein). Choosing for simphcity pp ~ y/3 
one has 

= V2(±x±y±z), (13) 

e<"^-±2xTyTz, (14) 

e<"' = V3(±yTz). (15) 

All triads can be transformed to each other by rotations and/or reflection. So in this 
model one obtains four identical copies of right and left relativistic Weyl fermions. 
They may be considered as analogs of generations of Standard Model fermions, 
but with rig = 4. Different, but related mechanism for the origin of generations is 
suggested in 1601 . 



4 Exotic fermions 

In many systems (including condensed matter and relativistic quantum vacua), the 
Fermi points with elementary charges = ±1 may merge together forming either 
the neutral point with = or point with multiple (i.e. |A^| > 1 ||6T| ). In this 
case topology and symmetry become equally important, because it is the symme- 
try which may stabilize the degenerate node. Example is provided by the Standard 
Model of particle physics, where 16 fermions of one generation have degenerate 
Dirac point at p = with the trivial total topological charge A^ = 8 — 8 = 0. In the 
symmetric phase of Standard Model the nodes in the spectrum survive due to a 
discrete symmetry between the fermions and they disappear in the non-symmetric 
phase forming the fully gapped vacuum ITI. In case of degenerate Fermi point with 
|A^| > 1, situation is more diverse. Depending on symmetry, the interaction between 
fermionic flavors may lead to splitting of the multiple Fermi point to elementary 
Dirac points ll56l : or gives rise to the essentially non-relativistic energy spectrum 
E± (/7 — > 0) — > ±p^, which corresponds to different scaling for space and time in the 
infrared: r — > ^r, f — > b^t. The particular case of anisotropic scaling with A^ = 3 was 
suggested by Hofava for quantum gravity at short distances, the so-called Hoi'ava- 
Lifshitz gravity ll62l 16311641 . while the anisotropic scaling in the infrared was sug- 
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gested in Ref. ||65l . The topology of the multiple Dirac point provides another pos- 
sible realization of anisotropic gravity, which is different from the scenario based 
on Lifshitz point in the theory of phase transitions ||66l . 



4.1 Dirac fermions with quadratic spectrum 

The non-linear spectrum arising near the Fermi point with N = 2 has been dis- 
cussed for different systems including graphene, double cuprate layer in high- 
Tc superconductors, surface states of topological insulators and neutrino physics 
II§7l[Tl|2l|29l|68l|69l|70l|7l]|71. The spectrum of (quasi)particles in the vicin- 
ity of the doubly degenerate node depends on symmetry. Let us consider zeroes of 
co-dimension 2 - point nodes in 2+1 systems. The node with topological charge 
= +2 takes place in bilayered graphene. According to general classification ifTSll . 
the topology alone cannot protect zeroes of co-dimension 2: the topological invari- 
ant takes place only in the presence of a symmetry. In particular, if we restrict con- 
sideration only to real (Majorana) fermions, the nodes obey Z2 topology with sum- 
mation law 1 + 1=0 ifTSl . To make the multiple Fermi point possible we need an 
additional symmetry K which extends the group Z2 to the full group of integers Z. 
The relevant symmetry protected topological invariant is 1731 |2l |29l l74l : 

N = tr i dl KG(o)=0,p)diG-\co = Q,p)=tr i dlKJ^-\p)di.y/f{p), 
Am JC JC 

(16) 

where C is contour around the Dirac point in 2D momentum space {p.x,Py)', K is 
the relevant symmetry operator; G is the Green's function matrix at zero frequency, 
which can be used as the effective Hamiltonian, ^(p) = G^' (o = 0,p); the oper- 
ator K commutes or anticommutes with the effective Hamiltonian. 

Provided the symmetry K is preserved and thus the summation law for takes 
place, one finds several scenarios of the behavior of the system with the total topo- 
logical charge = +2. 

(i) One may have two fermions with the linear Dirac spectrum, with the nodes 
being at the same point of momentum space. This occurs if there is some special 
symmetry, such as the fundamental Lorentz invariance. 

(ii) Exotic massless fermions emerge. In the 2D systems, these are gapless Dirac 
fermions with parabolic energy spectrum, which emerge at low energy: 

£±(p)«V- (17) 
They are described by the following effective Hamiltonian 

^bx,A.)= ((^^_°^^)2 + =(p.?-p?)cT,-2p,p,(T2, (18) 
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(a) (b) 



Fig. 6 Illustration of summation rule for momentum-space topological invariant. Splitting ofN = 2 
point with quadratic dispersion (a) into four Dirac points: A' = 2= l + 1 + 1— 1 (b). 



where (7i and (72 are Pauli matrices. The topological charge of the node at the 
point px = Py = is given by the Eq.(fT6]l. where the symmetry operator K is rep- 
resented by the Pauli matrix 03. With effective Hamiltonian ( fTSI ) one obtains that 
the node with the quadratic spectrum has the charge = 2. In 3D systems, the 
corresponding fermions with the topological charge = 2 in Eq. (|2]i are the semi- 
Dirac fermions, with linear dispersion in one direction and quadratic dispersion in 
the other Ig?) : 

£±(p)«±ycV?+pl- (19) 

(iii) The Weyl point with = +2 may split either into two Weyl points each with 
A^ = +1 (see 1561 for the relativistic 3+1 system) or into four Weyl points (three with 
N = +1 and one with N = — I, see Fig.|6]l. The effective Hamiltonian for the latter 
case is l|75ll76ll : 

The energy spectrum of this Hamiltonian has four Dirac point: the node at p = 
has the topological charge N = — I, while three nodes at px + ipy = —se^'^'^'l^ with 
integer k have charges A^ = + 1 each, so that the summation rule A^=:l + 1 + 1 — 1=2 
does hold. 

Note that in options (ii) and (iii) the (effective) Lorentz invariance of the Dirac 
point is violated. This suggests that the topological mechanism of splitting of the 
Dirac point ||56l or of the formation of the nonlinear dispersion ||67l may lead to 
the spontaneous breaking of Lorentz invariance in the relativistic quantum vacuum, 
which in principle may occur in the neutrino sector of the quantum vacuum 1,77. ,78l 

ED. 

Let us consider the Fermi point with higher degeneracy, described by the sym- 
metry protected topological invariant N > 2. 
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4.2 Dime fermions with cubic spectrum 

Let us consider first the case with N = 3 lISOl . Examples are three families of right- 
handed Weyl 2-component fermions in particle physics; three cuprate layers in high- 
Tc superconductors; three graphene layers, etc. If the Fermi point is topologically 
protected, i.e. there is a conserved topological invariant A^, the node in the spectrum 
cannot disappear even in the presence of interaction, but it can split into nodes 
with elementary charge N ~ I. The splitting can be prevented if there is a symmetry 
in play, such as rotational symmetry. Here we provide an example of such symmetry. 




-2 -1 1 2 

9 

Fig. 7 Splitting of Fermi point with N = 3 into three Fermi points with W = 1 in the model 
discussed in iSOj . 



We consider 3 species (families or flavors) of fermions, each of them being de- 
scribed by the invariant = +1 in Eq.(fT6]l and an effective relativistic Hamiltonian 
emerging in the vicinity of the Fermi point 

(P) = C7 • P = a,p_, + GyPy . (21) 

The matrix K = a, anticommutes with the Hamiltonian. This supports the topolog- 
ically protected node in spectrum, which is robust to interactions. The position of 
the node here is chosen at p = 0: 



(22) 
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The total topological charge of three nodes at p = of three fermionic species is 
= +3. Let us now introduce matrix elements which mix the fermions. If these 
elements violate symmetry K, then the topological invariant cannot be constructed 
and point node will be destroyed, so let the matrix elements obey the symmetry K. 
For the general case of the matrix elements, but still obeying the symmetry K, the 
multiple node will split into 3 or more elementary nodes, obeying the summation 
rule: 3 = 1 + 1 + 1 = 1 + 1 + 1 + 1 — 1 = - (see Fig.|7]l. However, in the presence of 
some extra symmetry, which prevents splitting, the branch with the cubic spectrum 
emerges. 



E/g 




























Fig. 8 Spectrum of the Hamiltonian <23b showing cubic dispersion for the lowest two eigenvalues 
around the point p = 0. The spectrum has been calculated with equal coupling strengths gi2 = 
Si3 ~ S2i = g, but cubic spectrum characterized by topological charge N = 3 preserves for any 
nonzero values of coupling strengths. Spectrum is shown as function of p, at p,. = 0. 

Example is provided by the following Hamiltonian ISOl 

^(P)= U21CJ- CJ-p g23C7+ , (23) 

where a* = ^ ((7,- ± iOy) are ladder operators. This Hamiltonian anti-commutes with 
K = G- and thus mixing preserves the topological charge in ( fTSI ). At = py = 
it is independent of the spin rotations up to a global phase of the coupling constants. 
Under spin rotation by angle Q all elements in the upper triangular matrix are mul- 
tiplied by e'^, while all elements in the lower triangular matrix are multiplied by 
e^'^ . This symmetry of triangular matrices does not allow the multiple Fermi point 
to split at /7 = 0, as a result the gapless branch of spectrum in Fig. |8]has the cubic 
form at low energy, £ -> 0, which corresponds to the topological charge A' = +3: 



Topology of quantum vacuum 



21 



In the low-energy limit the spectrum in the vicinity of the multiple Fermi point 
([24t is symmetric under rotations. But in general the spectrum is not symmetric 
as demonstrated in Figs. |8] There is only the symmetry with respect to reflection, 
{px,Py) —J- {Px,—Py)- The rotational symmetry of spectrum (l24li is an emergent phe- 
nomenon, which takes place only in the limit p 0. For trilayer graphene this spec- 
trum has been discussed in Ref. ISTI . 



4.3 Dirac fermions with quartic and higher order spectrum 



3 




2 












-2 




-3 





Fig. 9 Spectrum of the Hamiltonian ([25) showing quartic dispersion for the lowest two eigen- 
values around the point p = 0. Different colors correspond to different eigenvalues, spectmm is 
shown as function of p, at p,, = 0. The spectra have been calculated with equal coupling strengths 

gn = gl3 = g22 = gl4 = g24 = gu = g- 



In case of four fermionic species, the mixing which does not produce splitting 
of the Fermi point is obtained by the same principle as in (l23l i: all matrix elements 
aboye the main diagonal contain only (7+ (or (7^): 



^(p) 



/ c^ P 

821 <y^ 



8i2(y 

op 

832<y' 
842<y' 



8n<y 8i4(y^ 

823<y^ 824<y'^ 
842(y' (y 1? 



\ 



(25) 



Then again under spin rotation by angle 9 all elements in the upper triangular matrix 
are multiplied by e'^, while all elements in the lower triangular matrix are multiplied 
by e^'^ . As a result the multiple Fermi point with = 4 is preserved giving rise to 
the quartic spectrum in vicinity of the Fermi point: 
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For the tetralayer graphene this spectrum was suggested in Ref. ||82l . Again the rota- 
tional symmetry emerges new the Fermi point, but spectrum is not symmetric under 
permutations. In fact, there is no such 4x4 matrix that would consist of couplings 
described by ladder operators and which would be symmetric under permutations. 
In general the Fermi point with arbitrary may give rise to the spectrum 

£- = l^p'". (27) 

Such spectrum emerges in multilayered graphene 1291 l83l . The discussed symmetry 
of matrix elements g„,„ extended to 2N x 2N matrix gives ( |27] l with the prefactor 



\g\2\\g2-i\---\gN-\M 

Violation of this symmetry may lead to splitting of the multiple Fermi point into 
elementary Fermi points - Dirac points with A' = 1 and 'relativistic' spectrum 

£■2 oc p2. 

The effective Hamiltonian describing fermions in the vicinity of multiple Fermi 
point is H = + O^pl, see 1291 |3T1 . An example of the effective Hamiltonian 

describing the multiple Fermi point with topological charge in 3+1 systems is 

ED 

H = a,p, + a-pl + a+p^_ . (29) 

This Hamiltonian has the spectrum E'^ ~ pi + p^, which has linear dispersion in 
one direction and non-linear dispersion in the others. 

This is also apphcable to the vacuum of particle physics. The Lorentz symmetry 
prohibits both the splitting of the Dirac points and the non-linear non-relativistic 
spectrum. Situation changes if the Lorentz symmetry is viewed as an emergent phe- 
nomenon, which arises near the Dirac point (Fermi points with A^ = ± 1 ). In this case 
both both splitting of Dirac points and formation of non-linear non-relativistic spec- 
trum in the vicinity of the multiple Fermi point are possible, an the choice depends 
on symmetry. In both cases the mixing of fermions violates the effective Lorentz 
symmetry in the low-energy corner. This phenomenon, called the reentrant viola- 
tion of special relativity i67J , has been discussed for Nf ^ 3 fermion families in 
relation to neutrino oscillations ll58l . Influence of possible discrete flavor symme- 
tries on neutrino mixing has been reviewed in Ref. l85l . 



5 Flat bands in topological media 

Let us extend the results of the Sec. |4.3| to the limit case of infinite number of lay- 
ers, N ^ °° 1841 . If the layers are equivalent and interact only via nearest neighbor 
couplings, i.e. the nonzero matrix elements are 



gn — g23 — ■ --gN-lM = t , 



(30) 
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Px 




c 



N = 1 



Fig. 10 Nodal spiral generates topologically protecetd flat band on the surface. Projection of spiral 
on the surface determines boundary of flat band. At each (p., ,Pv) except the boundary of circle 
Px + Py = the system represents the ID gapped state (insulator). At each (px,Py) inside the 
circle, the insulator is topological being described by non-zero topological invariant and thus 
one has gapless edge state. The manifold of these zero-energe edge state inside the circle fonns the 
flat band found in Ea.l l32t . 



the low-energy spectrum becomes 



In the N ^ °° Umit one obtains that in the lowest energy branch all the fermions 
within circumference |p| = f have exactly zero energy. 



Here we consider the topological origin of the flat band. When the number of layers 
N ^ °°, the quasi two dimensional system transforms to the D = 3 system. We find 
that the flat band emerges in the first and the last layers, i.e. on the boundaries of 
the 3D system. Its formation is accompanied by the simultaneous formation of the 
line of node (zeroes of co-dimension 2) in the bulk material. These two objects, 
surface flat band and line of zeroes in bulk, are connected via a special kind of the 
topological bulk-surface correspondence: projection of the nodal line to the surface 
determines the boundary of the flat band (Fig.fTOli. 




(31) 



£(A^^oc,|p| <f) =0. 



(32) 
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5.1 Topological origin of surface fiat band 

To understand the topological origin of this branch and its structure let us consider 
the spectrum in the continuous limit. The effective Hamiltonian in the 3-dimensional 
bulk system which emerges in the limit of infinite number of layers is the following 
2x2 matrix 

Here t is the magnitude of the hopping matrix element between the layers in Eq.dSOl). 
The energy spectrum of the bulk system 

E'^ = [px-fcos(fl/5-)]^ + bv + fsin(fl/?,)]2, (34) 

has zeroes on line (see Fig.fTOli: 

Px = f cos(ap-) , py = —tsm{ap:) , (35) 

which forms a spiral, the projection of this spiral on the plane p- = const being the 
circle = p^ + Py = t^. This nodal line is topologically protected by the same 
topological invariant as in Eq.(fT6]l 

N ^ — tr <£ dl a-H-^ViH, (36) 
47ti Jc 

where the integral is now along the loop C around the nodal line in momentum 
space, see Fig. [TO] The winding number around the element of the nodal line is 

N=l. 

Let us consider now the momentum px as a parameter of the ID system, then 
for IpxI 7^ t the system represents the fully gapped system - ID insulators. This 
insulator can be described by the same invariant as in Eq.([36]l with the contour of 
integration chosen parallel to p-, i.e. along the ID Brillouin zone at fixed px (due to 
periodic boundary conditions, the points p^ ^±n/a are equivalent and the contour 
of integrations forms the closed loop): 

A^(P±) = r-t'" / dp,a,H-Wp_H. (37) 

4-711 J-K/a 

For |px| < t the ID insulator is topological, since A^(|pxl <t) = 1, while for |px| > t 
one has A^i (px ) = and the ID insulator is the trivial band insulator The line |px | = 
/ thus marks the topological quantum phase transition between the topological and 
non-topological ID insulators. 

Topological invariant A^(px) in ( |37] ) determines the property of the surface bound 
states of the ID system at each px. Due to the bulk-edge correspondence, the topo- 
logical ID insulator must have the surface state with exactly zero energy. Since such 
states exist for any parameter within the circle |px | ~ t, one obtains the flat band of 
surface states with exactly zero energy, £ ( | Px | <t)=0, which is protected by topol- 
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Fig. 11 Formation of the surface flat band. When the number A' of layers increases, the disper- 
sionless band evolves from the gapless branch of the spectrum, which has the form E = ±|p^|'^ 
in the vicinity of multiple Dirac point. The spectrum is shown as a function of for Py = 0. The 
curves for A' = 100 and A' = 200 are almost on top of each other. Asymptotically the spectrum 
E = ±|pj^|'^ transforms to the dispersionless band within the projection of the nodal line to the 
surface. 

ogy. This is the origin of the unusual branch of spectrum in Eq.(l32b: it represents the 
band of topologically protected surface states with exactly zero energy. Such states 
do not exist for parameters |px | > t, for which the ID insulator is non-topological. 

The zero energy bound states on the surface of the system can be obtained di- 
rectly from the Hamiltonian: 

H = Ox{px-tco^{ap,)) + ay{py + t^\n{ap;)) , p^ = —id, , z<0. (38) 

We assumed that the system occupies the half-space z < with the boundary at 
Z = 0. This Hamiltonian has the bound state with exactly zero energy, ^(px) = 0, 
for any |px| < t, with the eigenfunction concentrated near the surface: 



The normalizable wave functions with zero energy exist only for px within the circle 
IPxI < U i-C- the surface flat band is bounded by the projection of the nodal spiral 
onto the surface. Such correspondence between the flat band on the surface and lines 
of zeroes in the bulk has been also found in Ref. 1321 for superconductors without 
inversion symmetry. 




(39) 
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Fig. 12 Fomiation of the nodal line from the evolution of the gapped branch of the spectrum of the 
multilayered system, when the number N of layers increases. The spectrum is shown as function of 
Px for py = 0. The curves for N = 100 and N = 200 lie almost on top of each other, indicating the 
bulk limit. Asymptotically the nodal line p, = fcos(ap.), Py = — fsin(ap;) is formed (two points 
on this line are shown, which correspond to Py = 0). 



5.2 Dimensional crossover in topological matter: Formation of the 
flat band in multilayered system 



The discrete model with finite number of layers has been considered in Ref . Il84ll . 
It is described by the 2N x 2N Hamiltonian with the nearest neighbor interaction 
between the layers in the form: 

////px) = (T-p±5,7-f(7+5,j+i -fc7"5,j_i , l<i<N , Px = (/7^,/7v). (40) 

In the continuous limit of infinite number of layers ( l40b transforms to (l33T l with 
the nodal line in the spectrum, while for finite the spectrum contains the Dirac 
point with multiple topological charge equal to A^. Figures (fTTT i - (fl4] | demonstrate 
how this crossover from 2D to 3D occurs. When the number A^ of layers increases, 
the dispersionless surface band evolves from the gapless branch of the spectrum 
E = ijpxl'^. Simultaneously, the gapped branches of the spectrum of the finite-A^ 
system give rise to the nodal line in bulk. This scenario of formation of the surface 
flat band takes place if the symmetry does not allow the splitting of the multiple 
Dirac point, and it continuously evolves to the dispersionless spectrum. 

If this symmetry is absent, but the symmetry supporting the topological charge 
persists, the scenario of the flat band formation is different but still is governed by 
topology. Fig. ( fTsT i demonstrates the formation of the flat band in this situation. The 
reason, why the flat band emerges is the formation of the topologically protected 
nodal line in bulk. The projection of the nodal line on a surface gives the boundary 
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Fig. 13 The lowest energy states for different p., and p,. and arbitrary p- in the bulk limit. The flat 
band of surface states is formed in the region p\+p^, <t^. Outside this region only the bulk states 
exist. In a simple model considered here flat band comes from the degenerate Dirac point with 
nonlinear dispersion. However, this is not necessary condition: the flat band emerges whenever the 
nodal line appears in the bulk. 




Py/i 

Fig. 14 Result of the transformation of the gapped states in the process of dimensional crossover. 
They form the nodal line in bulk p, =?cos(op.), p,, = — f sin(ap-) whose projection to the (p, , p, )- 
plane is shown. 



of the flat band emerging on this surface. This is the realization of the bulk-surface 
correspondence in systems with the nodal lines in bulk. The other examples can be 
found in Ref. US). 
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Fig. 15 Evolution of the spectrum in a different model, in which the multiple Dirac point is absent, 
but the suiface flat band is formed together with the formation of the nodal line in bulk 

6 Anisotropic scaling and Hof ava gravity 

We know that in the vicinity of the Weyl point with elementary topological charge 
A^ = +l orA^ = — 1 the quantum electrodynamics and gravity emerge as effective 
fields. Here we discuss what kind quantum electrodynamics emerges near the node 
with higher A^. This is the quantum electrodynamics with anisotropic scaling, in 
which as distinct from isotropic relativistic QED the electric and magnetic fields 
obey different scaling laws. The analogous anisotropic scaling has been discussed 
for quantum gravity at short distances. 



6.1 Effective theory near the degenerate Dirac point 

Let us consider again the 2+1 system with the multiple Dirac point and write the 
following effective Hamiltonian: 



J^N = ((ei +/e2) • (p-eA))l"l + ((ei - /e2) • {v-eP^)P . 



For a single layer (A^ = 1) the Hamiltonian ( |4TI ) is reduced to the conventional Weyl- 
Dirac Hamiltonian for massless particles in 2+1 dimension; 



=C7-'ei-(p-eA) + (7-^'e2-(p-eA)=4c7"(/7;-eA,) , fl = (l,2). (42) 



Here the vectors ei and 62 play the role of zweibein (in the ground state they are 
mutually orthogonal). The vector A is either the vector potential of the effective 
electromagnetic field which comes from the shifts of the node, or the real electro- 



(41) 
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magnetic field, as it takes place for electrons in graphene. The whole dreibein e„ 
with a ~ (1,2,3) and = (0, 1,2) emerges for the Green's function and gives the 
effective metric 2+1 metric as a secondary object: 

g^' = il"'e^el , (43) 

For general |A^| situation is somewhat different. While the action for the relativis- 
tic fermions is invariant under rescaling r = br', t = bt' , the action for the fermions 
living in the vicinity of the multiple Dirac point is invariant under anisotropic rescal- 
ing r = br' , t = Zjl^lf'. The anisotropic scaling is in the basis of the Hofava gravity, 
which is described by the space components of metric are separated from the time 
component and have different scaling laws [62 63 MllSl- The square of Hamilto- 
nian (4T[ gives the space metric in terms of zweibein: 

jr^2^4=(g'^(p,-M,)(p,-M,))'"' , g'^ =e\ei+e^ei. (44) 



6.2 Effective electromagnetic action 

If the effective action for fields Ca and is obtained by integration over fermions in 
the vicinity of the multiple Dirac point, this bosonic action (actually the terms in ac- 
tion which mostly come from these fermions) inherits the corresponding conformal 
symmetry of the massless fermions. For anisotropic scaling the conformal invari- 
ance means invariance under g''^ — > b^g''' and g^^ — > b^^^^g^^; y/—g 
(where D is space dimension); while g^' is not considered. 

6.2.1 Single layered graphene and relativistic fields 

For D y^2 the spectrum of multiple Fermi point becomes more complicated, and in 
general is not isotropic, see Eq.(l29]l. For general D the spectrum is isotropic only 
for |A^| = 1, where one obtains effective relativistic massless D+l quantum electro- 
dynamics, which is Lorentz invariant. This implies the following nonlinear action 

Se,„{\N\ 1,£>) = j d'^xdt [B^ _£2j ^ ^ (-45-) 

For D = 3 the action is proportional to (Z?^ — E^)\n{B^ — E^), and is imaginary at 
< giving rise to Schwinger pair production in massless quantum electrody- 
namics. The similar imaginary action tales place for < E^ for D ^3. For exam- 
ple, for a single layer graphene (D = 2, |A^| = 1) reproduces the relativistic 2+1 QED 
which gives rise to Lagrangian (B^ — fi^)^/^ fST! with the running coupling constant 
1/a = {3/ 2)/ 871^. The action is imaginary at B^ < E^ which corresponds to 
Schwinger pair production with the rate E^^^ at B = 0. 
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6.2.2 Bilayer graphene 

For bilayered graphene, assuming the quadratic dispersion |A^| = D = 2, the expected 
conformal invariant Heisenberg-Euler action for the constant in space and time elec- 
tromagnetic field, which is obtained by the integration over the 2+ 1 fermions with 
quadratic dispersion, is the function of the scale invariant combination jj. Il88l[89l : 

£■2 

d^xdt B^gipL) , ^ = (46) 

The asymptotical behavior in two limit cases, ^(/i 0) ~ const and — > oo) ~ 
^2/3^ gives the effective actions for the constant in space and time magnetic and 
electric fields: 



Sb 



I J dhdt , SE^{b + ic) j d^xdtE^I^ . (47) 



The parameter a is the logarithmic coupling constant; the parameter b describes the 
vacuum electric polarization; and the parameter c describes the instability of the 
vacuum with respect to the Schwinger pair production in the electric field, which 
leads to the imaginary part of the action. The action also contains the linear non- 
local term 

Snon-locaim ^ D ^ 2) = J d" xdt ^g^' g'" Fak—^FQ„ , (48) 

which corresponds to the polarization operator 

7100--== , \N\=D = 2. (49) 
V^* - or- 



6.2.3 D=2 systems with nodes with topological charge A' 

In general case of a 2D system with A^-th order touching point in spectrum (a kind of 
multilayered graphene) the Heisenberg-Euler action contains among the other terms 
the following nonUnear terms in the actions for magnetic and electric fields 



Sb{N,D^2)^ Jd^xdtB^ , Se{N,D = 2)^ j d^xdt{-E^Y^(^+W) , (50) 

where imaginary part of the action is responsible for the pair production in electric 
field |i88„90J, while the Hnear action of the type ( |48] | corresponds to the polarization 
operator 
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The expected action for gravitational field is 
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Sgrilv — d xdt^/—g 



(52) 



where Ki,K2,... are dimensionless quantities, and the other terms of that type are 
implied. 



7 Fully gapped topological media 



Examples of the fully gapped topological media are the so-called topological band 
insulators in crystals [3]. Examples are Bi2Se3, Bi2Te3 and Sb2Te3 compounds 
which are predicted to be 3+1 topological insulators |91j. But the first discussion 
of the 3+1 topological insulators can be found in Refs. 1921 |93l . The main feature 
of such materials is that they are insulators in bulk, where electron spectrum has 
a gap, but there are 2+1 gapless edge states of electrons on the surface or at the 
interface between topologically different bulk states as discussed in Ref. ll93l . The 
similar properties are shared by the fully gapped 3D topological superfluids and su- 
perconductors. The spin triplet p-wave superfluid ^He-B represents the fully gapped 
superfluid with nontrivial topology. It has 2+1 gapless quasiparticles living at inter- 
faces between vacua with different values of the topological invariant describing 
the bulk states of ^He-B lll4ll94ll . The quantum vacuum of Standard Model below 
the electroweak transition, i.e. in its massive phase, also shares the properties of the 
topological insulators and gapped topological superfluids and is actually the rela- 
tivistic counterpart of ■'He-B 1231 . 

Examples of the 2+1 topological fully gapped systems are provided by the films 
of superfluid ^He-A with broken time reversal symmetry 12111951 and by the planar 
phase which is time reversal invariant l2ni95]| . The topological invariants for 2+1 
vacua give rise to quantization of the Hall and spin-Hall conducticity in these films 
in the absence of external magnetic field (the so-called intrinsic qauntum and spin- 
quantum Hall effects) 12111961 . 



7.1 2+1 fully gapped vacua 

7.1.1 ^He-A film: 2+1 chiral superfluid 

The gapped (nodeless) ground states (vacua) in 2+1 systems are characterized by 
the invariant obtained by dimensional reduction from the topological invariant de- 
scribing the nodes of co-dimension 3. The invariant for the Fermi point in (|2|is the 
integral over the 3D surface a around the singularity. The invariant describing the 
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Fig. 16 Skyrmion in p-space with momentum space topological charge A' = — 1 in (55). It de- 
scribes topologically non-trivial vacua in 24-1 systems with a fully gapped non-singular Green's 
function. Vacua with nonzero N have topologically protected gapless edge states. The nonzero 
topological charge leads also to quantization of Hall and spin Hall conductance. 



2D insulator or nodeless supefluid is the integral over the whole (2+l)-dimensional 
momentum-frequency space {px-Py, co): 



24n- 



d^pdco Gdp. G-^GdpG-^ Gdp^ G ' 



(53) 



This equation ( l53T l was introduced in relativistic 2-1-1 theories ifTTl [121 [131 and 
for the film of ^He-A in condensed matter |[T6l [2TI . where it was inspired by the 
dimensional reduction from the Fermi point, see 1951 . In simple case of the 2x2 
matrix, the Green's function can be expressed in terms of the three-dimensional 
vector d(/7^,/7y), 

G-\(0,px.p,)^i(0-H , H^x-d{p,„py), (54) 

where H is the Bogoliubov - de Gennes Hamiltonian for fermions in p-wave super- 
fluids, and T are the Pauli matrices. Example of the d-vector configuration, which 
corresponds to the topologically nontrivial vacuum is presented in Fig. [16] This is 
the momemtum-space analog of the topological object in real space - skyrmion. In 
real space, skyrmions are described by the relative homotopy groups 1971 ; they have 
been investigated in detail both theoretically and experimentally in the A phase of 
^He, see Sec. 16.2 in 11] and the review paper 1981 . 

For the Green's function in i54[ the winding number of the momentum-space 
skyrmion in Eq.([53]l is reduced to lT6l 

N^T- d^pd-i^x — ] , (55) 



47t J \dpx dpy 



where d = d/|d| is unit vector. For a single layer of the ^He-A fibn and for one spin 
projection, the simplified Bogoliubov - de Gennes Hamiltonian has the form: 
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Fig. 17 Dependence of the topological invariant i53\ on the thickness of ^^He-A film. The even 
values of N result from the spin degeneracy. At the topological phase transitions between the states 
with different A', the gap in the spectrum of fermions is nullified. The interface between the vacua 
with different topological charges contains gapless fermions - edge states. The number of the 
gapless fermions is related to the difference of the topological charges by index theorem. 



H = T- d(p) = iCO + T3 y-^^^ - M j + + TlPy , (56) 

For ;Li > the topological charge in Eq.dSSTi or in Eq.(l53]l is = 1, and the d(p) 
field forms the skyrmion in momentum space of the type shown in Fig. [16] For 
jj. < the topological charge is trivial, N = 0. That is why at /i = there is a 
topological quantum phase transition between the topological superfluid at /i > 
and non-topological superfluid at /i < 1951 . 

In general case of multilayered ^He-A, topological charge A^ may take any integer 
value of group Z. This charge determines quantization of Hall and spin-Hall conduc- 
tance and the quantum statistics of the topological objects - real-space skyrmions 
|[l6l|2ll|96l|95l. For A^ = 4;t + 1 and A^ = 4yt + 3, skyrmion is anyon; for A^ = 4yt + 2 
it is fermion; and for A^ = 4A: it is boson ||95|| . This demonstrates the importance 
of the Z2 and Z4 subgroups of the group Z in classification of topological matter; 
and also provides an example of the interplay of momentum-space and real-space 
topologies. 



7.1.2 Planar phase: time reversal invariant gapped vacuum 



In case when some symmetry is present, additional invariants appear, which corre- 
spond to dimensional reduction of topological invariant supported by symmetry Nk 
in®: 

-1 



j cf-pd(0 KGdpfi- ' Gdpj G ' Gdp^G- 



(57) 
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where the matrix K commutes with the Green's function matrix. Example of the 
symmetric 2+1 gapped state with Nk is the film of the planar phase of superfluid 
^He 1211 |95l . In the single layer case, the simplest expression for the Green's func- 
tion is 



where H is the Bogoliubov - de Gennes Hamiltonian for fermions in this spin-triplet 
p-wave superfluid. The symmetry operator K, which supports the topological in- 
variant \sK ~ 13(7-: it commutes with the Green's function. The planar state is time 
reversal invariant. It has trivial conventional topological charge N ^0 and non-zero 
symmetry protected charge Nk ~ 2. For the general case of the quasi 2D film with 
multiple layers of the planar phase, the invariant Nk belongs to the group Z. The 
magnetic solid state analog of the planar phase is the 2D time reversal invariant 
topological insulator, which experiences the quantum spin Hall effect without exter- 
nal magnetic field ||3l . 



8 Relativistic quantum vacuum and superfluid ^He-B 

Let us now turn to the class of 3+1 fully gapped systems, which is represented by 
Standard Model in its massive phase and superfluid ^He-B. 

In the broken symmetry phase of Standard Model below the electroweak transi- 
tion, symmetry does not support the topological invariant responsible for the nodes 
in spectrum. In this phase there is no mass protection by topology and thus all the 
fermions become massive, i.e. Standard Model vacuum becomes the fully gapped 
insulator 

In quantum liquids, the fully gapped three-dimensional system, which is simi- 
lar to the vacuum of massive Standard Model, is represented by another phase of 
superfluid ■'He - the ^He-B. This phase has time reversal symmetry and nontrivial 
topology supported by symmetry, which gives rise to the 2D gapless quasiparticles 
living at interfaces between vacua with different values of the topological invariant 
or on the surface of ^He-B |[T4l |94l |99l [TOOl . 



8.1 Superfluid ^He-B 

^He-B belongs to the same topological class as the vacuum of Standard Model in 
its present insulating phase ||23| . The topological classes of the ^He-B states can be 
represented by the following simplified Green's function and the Bogoliubov - de 
Gennes Hamiltonian: 




(58) 
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G (co,p) = i(o-H , // = T3 M I + Tic^a p. 

' 2m 



(59) 



In the limit 1/m = this model ^He-B transforms to the vacuum of massive rela- 
tivistic Dirac particles with speed of light c = and mass parameter M = —ji. 

In the fully gapped systems, the Green's function has no singularities in the whole 
4-dimensional space (ft),p). That is why we are able to use the Green's function 
at (B = 0, which corresponds to the effective Hamiltonian, i/eff(p) = ^G^'(0,p). 
The topological invariant relevant for ^He-B and for quantum vacuum with massive 
Dirac fermions is: 



247r2 



tr 



(B=0 



d^'p K Gdpfi- ' Gdp G- ' Gdpfi 



— 1 



(60) 



with matrix = T2 which anti-commutes with the Green's function at = 0. In 
^He-B, the T2 symmetry is combination of time reversal and particle-hole symme- 
tries; for Standard Model the matrix T2 = 757". Note that at CO = the symmetry 
of the Green's function is enhanced, and thus there are more matrices K, which 
commute or anti-commute with the Green's function, than at O) 7^ 0. 
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Fig. 18 Phase diagram of topological states of ■'He-B in Eq. <59| l in the plane [p., 1/m). States on 
the line 1 /m = correspond to the Dirac vacua, which Hamiltonian is non-compact. Topological 
charge of the Dirac fermions is intemiediate between charges of compact ^He-B states. The line 
l/m — separates the states with different asymptotic behavior of the Green's function at infinity: 
G^' ((B = 0,p) — > ±tifp- /2m. The line = marks topological quantum phase transition, which 
occurs between the weak coupling ■'He-B (with > 0, ;;j > and topological charge Nk = 2) 
and the strong coupling ■'He-B (with /i < 0, m > and Nk = 0). This transition is topologically 
equivalent to quantum phase transition between Dirac vacua with opposite mass pai'ameter M = 
it|/i|, which occurs when /i crosses zero along the line 1/m = 0. The interface which separates 
two states contains single Majorana fermion in case of ^He-B, and single chiral fermion in case 
of relativistic quantum fields. Difference in the nature of the fermions is that in Fermi superfluids 
and in superconductors the components of the Bogoliubov-Nambu spinor ai'e related by complex 
conjugation. This reduces the number of degrees of freedom compared to Dirac case. 



36 



G.E. Volovik 



Fig. [18] shows the phase diagram of topological states of ^He-B in the plane 
{jj., 1 /m). The line 1/m = corresponds to the Dirac vacuum of massive fermions, 
whose topological charge is determined by the sign of mass parameter M = —jj.: 

Nk = sign{M). (61) 

The real superfluid ^He-B lives in the corner of the phase diagram jj.>0, m > 0, 
fi ^ mc\, which also corresponds to the limit of the weakly interacting gas of ^He 
atoms, where the superfluid state is described by Bardeen-Cooper- Schrieffer (BCS) 
theory. However, in the ultracold Fermi gases with triplet pairing the strong coupling 
limit is possible near the Feshbach resonance 11011 . When pL crosses zero the topo- 
logical quantum phase transition occurs, at which the topological charge AOs: changes 
from Nk = 2 to Nk = 0. The latter regime with trivial topology also includes the 
Bose-Eistein condensate (BEC) of two-atomic molecules. In other words, the BCS- 
BEC crossover in this system is always accompanied by the topological quantum 
phase transition, at which the topological invariant changes. 

There is an important difference between ^^He-B and Dirac vacuum. The space of 
the Green's function of free Dirac fermions is non-compact: G has different asymp- 
totes at IpI oo for different directions of momentum p. As a result, the topologi- 
cal charge of the interacting Dirac fermions depends on the regularization at large 
momentum. ^He-B can serve as regularization of the Dirac vacuum, which can be 
made in the Lorentz invariant way 1231 . One can see from Fig. [18] that the topo- 
logical charge of free Dirac vacuum has intermediate value between the charges of 
the ^He-B vacua with compact Green's function. On the marginal behavior of free 
Dirac fermions see Refs. |[T5] [19] [111991 . 

The vertical axis separates the states with the same asymptote of the Green's 
function at infinity. The abrupt change of the topological charge across the line, 
ANk = 2, with fixed asymptote shows that one cannot cross the transition line adi- 
abatically. This means that all the intermediate states on the line of this QPT are 
necessarily gapless. For the intermediate state between the free Dirac vacua with 
opposite mass parameter M this is well known. But this is applicable to the gen- 
eral case with or without relativistic invariance: the gaplessness is protected by the 
difference of topological invariants on two sides of transition. 



8.2 From superfluid relativistic medium to ^He-B 

We have already seen, that the Dirac vacuum of massive Standard Model particles 
has nontrivial topology. As a result the domain wall separating vacua with opposite 
signs of the mass parameter M contains fermion zero modes II102I . Other examples 
of the topologically nontrivial states in relativistic theories can be provided by dense 
quark matter, where chiral and color superconductivity is possible. The topological 
properties of such fermionic systems have been recently discussed in Ref. 110311 . In 
particular, in some range of parameters the isotropic triplet relativistic superconduc- 
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tor is topological and has the fermion zero modes both at the boundary and in the 
vortex core. On the other hand, there is a range of parameters, where this triplet su- 
perconductor is reduced to the non-relativistic superfluid ^He-B 11041 . That is why 
the analysis in Ref. 11031 is applicable to ^He-B and becomes particularly useful 
when the fermions living in the vortex core are discussed. 

In relativistic superconductor or superfluid with the isotropic pairing - such as 
color superconductor in quark matter - the fermionic spectrum is determined by 
Hamiltonian 

H = T3{ca-p + pM- iir) + TiA , (62) 
for spin singlet pairing, and by Hamiltonian 

// = T3(ca-p + j3M-jLiff) + 75Ti4, (63) 

for spin triplet pairing II1041 11031 . Here a', jS and 75 are Dirac matrices, which in 
standard representation are 

M is the rest energy of fermions; is their relativistic chemical potential as dis- 
tinct from the non-relativistic chemical potential jj.; T„ are matrices in Bogoliubov- 
Nambu space; and A is the gap parameter. 

In non-relativistic limit the low-energy Hamiltonian is obtained by standard pro- 
cedure, see e.g. 11051 . The non-relativistic limit is determined by the conditions 

cp < M (65) 

and 



\M-^J^l + A^\<^M. (66) 

Under these conditions the Hamiltonian ( l62b reduces to the Bogoliubov - de Gennes 
(BdG) Hamiltonian for fermions in spin-singlet i-wave superconductors, while ( l63l l 
transforms to the BdG Hamiltonian relevant for fermions in isotropic spin-triplet 
/7-wave superfluid ^He-B in Eq. 



, ' \ R M „ A 



where the nonrelativistic chemical potential jj. = + A^ — M. 
The Dirac-BdG system in Eq.(l63]l has the following spectrum 



e = ±]jM^+c^p^ + A^ + lil±2^M^{nl + A-)+n^c^-p^. (68) 

This spectrum is plotted in Fig(T9] Depending on the value of the parameters Hr, A, 
M the spectral branches have different configurations. 



38 



G.E. Volovik 




Fig. 19 Plot of the spectrum of relativistic Hamiltonian <63b for two generic cases: (a) \i1r\ > jl^ 
when the minimum in the energy spectrum is away from the origin and (b) < fl^ when the 
minimum in the energy spectrum is at p = 0. At > /i^ there is a soft quantum phase transition 
between these two vacua. This transition is not topological, and thus the gap in the energy spectrum 
does not close at the transition. The gap closes at the topological transition occurring at fi^ = 
— as shown in plot (c). 



There is a soft quantum phase transition, at which the position of the minimum of 
energy E(p) shifts from the origin p = 0, and the energy profile forms the Mexican 
hat in momentum space. This momentum-space analog of the Higgs transition ||2l 
occurs when the relativistic chemical potential exceeds the critical value 

^*^=[^ + \l^+^^') ' ■ ^^^^ 

Figures [19] (a) and[T9](b) demonstrate two generic cases: l/i^j > /i^ when there are 
extremums of function e(p) at /9 7^ and \plr\ < jJ-g when all extremums are at 
the point p = 0. The formation of the Mexican hat at = /i^ is an example of 
non-topological quantum phase transition, which occurs without change of the topo- 
logical invariants and thus is not accompanied by the gapless intermediate state. Let 
us turn to the topological quantum phase transitions, at which the p-space topolog- 
ical invariant changes and the gap closes at the transition point as is shown in Fig. 
[11(c). 



8.3 Topology of relativistic medium and ^He-B 

Fig. [20] shows the phase diagram of the vacuum states of relativistic triplet super- 
conductors. Different vacuum states are characterized by different values of the 
topological invariant Nk in Eq.(i60]l, where the Green's function matrix at zero fre- 
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quency G^' (o) = 0,p) is equivalent to effective Hamiltonian. For ^He-B in Eq.(l67]i 
and for triplet relativistic superconductor in Eq. (|63] | the relevant matrix = T2, 
which anti-commutes with the Hamiltonian. The vacuum states with different Nk 
cannot be adiabatically connected, and thus at the phase transition lines the states 
are gapless. The circle /xj + 4^ = is an example of the line of topological quan- 
tum phase transition. In non-relativistic limit this corresponds to the line /i = in 
Figim The states inside the circle + = are topologically trivial, while the 
states outside this circle represent topological superconductivity II103II . The vacuum 
states with /ij + 4^ > and /ij +4^ < cannot be adiabatically connected 
which leads to the gap closing in Fig[T9](c). Discontinuity in the topological charge 
across the transition induces discontinuity in the energy of the ground state across 
the transition. For example, for the 2+1 px + ipy superfluid/superconductor the quan- 
tum phase transition is of third order, meaning that the third-order derivative of the 
ground state energy is discontinuous II106I . 




Fig. 20 Phase diagram of ground states of relativistic triplet superfluid in Eg. <63t in the plane 
{Hr,A). Topological quantum phase transitions are marked by thick lines. The states inside the 
circle pL^ + A- = M- are topologically trivial. The states outside this circle represent topological 
superconductors. The states on the lines of topological quantum phase transition are gapless. 



9 Fermions in the core of strings in topological materials 

In relativistic theories there is an index theorem which relates the number of fermion 
zero modes localized on a vortex with the vortex winding number |10?1. We know 
that the Dirac vacuum considered in Ref. 11071 has nonzero topological charge. 
This suggests that the existence of zero energy states in the core is sensitive not only 
to the real-space topological charge of a vortex, but also to the momentum-space 
topological charge of the quantum vacuum in which the vortex exists, and if so the 
index theorem should be extended to vortices in any fully gapped systems, including 
the non-relativistic superfluid ^He-B. Here we discuss this issue of the connection 
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between the topological charge of the vacuum and existence of the fermion zero 
modes on the topological objects in this vacuum. Our examples demonstrate that 
during the topological quantum phase transition from the topologically nontrivial 
vacuum to the trivial one, the fermion zero modes on vortices (or other objects) 
disappear. 



9.1 Vortices in ^He-B and relativistic strings 

Since in this problem both the momentum-space topology of bulk state and the 
real-space topology of the vortex or other topological defects are involved, the com- 
bined topology of the Green's function in the coordinate-momentum space (©, p, r) 
ll38l 11081 [T] 11091 II 101 seems to be relevant. However, though the bulk- vortex cor- 
respondence does evidently exist, the explicit index theorem which relates the exis- 
tence of the fermion zero modes to the topological charge of the bulk state and the 
vortex winding number is still missing. The existing index theorems are applicable 
only to particular cases, see e.g. II1081 11031 II 101 [TTTI . There is also a special in- 
dex theorem for superconductors/superfluids with a small gap A pL. Spectrum of 
fermions in these superconductors has branches which cross zero energy as a func- 
tion discrete quantum number - angular momentum L 11 121 . The index theorem 
relates the number of such branches with the vortex winding number 11 131 . Here 
we are interested in the true fermion zero modes - the branches of spectrum £(/?-), 
which cross zero as function of momentum along the vortex line. 

Example, which demonstrates that the connection between the topological charge 
Nk and the existence of Majorana fermions - fermion zero modes on vortices - is 
in Fig. |2T] For ^He-B, which lives in the BCS range of parameters where Nk ^ 0, 
the gapless fermions in the core have been found in Ref. [ 114|. On the other hand, 
in the strong coupling limit the ^He-B transforms to the Bose-Einistein condensate 
(BEC) of molecules. The latter does not contain fermionic excitations, and thus one 
should not expect the existence of gapless fermions in the vortex core. Thus one ex- 
pects that somewhere in the region of crossover between the BCS-like and BEC-like 
regime the spectrum of fermions localized on vortices must be reconstructed. On the 
other hand, the topological reconstruction of the fermionic spectrum in the vortex 
core cannot occur during the adiabatic deformation. The discontinuous deformation 
of the spectrum in the core is only possible during the topological quantum phase 
transition in bulk: at such transition the intermediate bulk gapless state is crossed 
which destroys the adiabaticity. This is just what happens: in the BEC limit, the 
chemical potential /i is negative and the topological charge of the vacuum state is 
trivial, Nk = 0. The reconstruction of the bulk spectrum at the topological quan- 
tum phase transition occuring at /i = 0, see Fig. [18] triggers reconstruction of the 
spectrum of fermion zero modes in the core. At < the topological charge Nk 
nullifies and simultaneously the gap in the spectrum of core fermions arises, see 
Fig. |2T| The similar situation, when the quantum phase transition in bulk produces 
leads to appearance or disappearance of fermion zero modes, has been discussed in 




Fig. 21 Schematic illustration of spectrum of the fermionic bound states in the core of the 
most symmetric vortex with « = 1, the so-called o- vortex 1981 . in fully gapped spin triplet su- 
perfluid/superconductor of ■'He-B type, (left): Spectrum of bound state in the ■'Hc-B, vortex which 
corresponds to the weak coupling BCS regime with non-zero topological charge Nk = 2 11141 . L is 
the azimuthal quantum number of fermions in the vortex core. There are two fermion zero modes, 
which cross zero energy in the opposite directions, (right): The same vortex but in the topologically 
trivial BEC regime, = 0, does not have fermion zero modes. The spectrum of bound states is 
fully gapped. Fermion zero modes disappear at the topological quantum phase transition, which 
occurs in bulk liquid at ;U = 0. Similar situation may take place for strings in color superconductors 
in quark matter 11031 . 



Ref. Ill 151 for the other type of p-wave vortices; in Ref. II103I for strings in color 
superconductors in quark matter; and in Refs. 01 16111 17l for Majorana fermions on 
the edges of quantum wire (the review on Majorana fermions in superconductors 
can be found in Ref. 11 181 ). 

Another example is provided by the fermions on relativistic vortices in Dirac 
vacuum discussed in Ref. II107II . The Dirac vacuum has the nonzero topological in- 
variant, Nk = ±1, see Fig. [18] This is consistent with the existence of the fermion 
zero modes on vortices, found in Ref. 11071 . The index theorem for fermion zero 
modes on these vortices can be derived using the topology in combined coordi- 
nate and momentum space. The number of fermion zero modes on a vortex A'zm 
can be expressed via the 5-form topological invariant in terms of Green's function 
G(£0,p,r) EOlinSl 
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(70) 

The space integral is along the closed contour C around the vortex line. For the 
vortex in Dirac vacuum, equation dTOb reproduces the index theorem discussed in 
Ref. II107II : the algebraic number of fermion zero modes equals the vortex winding 
number A^zm = «• 

For vortex in ^He-B one obtains Nzm ~ 0. This is, however, consistent with 
Fig. |2T] two branches of zero modes have opposite signs of velocity — dE / d p,. 
Though the algebraic sum of zero modes, A^zm = 1 — 1=0, due to the special sym- 
metry of the vortex configuration the two branches of fermion zero modes do not 
cancel each and the gap is not formed. To resolve the fermion zero modes in the sys- 
tems with symmetry, the index theorem for the zero modes must be supplemented 
by symmetry consideration. 

The 5-form topological invariant similar to Eq.dTOli has been discussed also in 
Ifn il20l . In particular it is responsible for the topological stability of the 3h-1 chiral 
fermions emerging in the core of the domain wall separating topologically differ- 
ent vacua in A+\ systems (see Sec. 22.2.4 in yj). The topological invariant for 
the general 2n + 1 insulating relativistic vacua and the bound chiral fermion zero 
modes emerging there have been considered in 01211 11221 l60l . Application of the 
5-form topological invariant to the states in lattice chromodynamics can be found in 



9.2 Flat band in a vortex core: analog of Dirac string terminating 
on monopole 



E(Pz) 



continuous spectrum 




Pz 



Fig. 22 Schematic illustration of the spectrum of bound states E(p~) in the vortex core. The 
branches of bound states terminate at points where their spectrum merges with the continuous 
spectrum in the bulk. The flat band terminates at points where the spectrum has zeroes in the bulk, 
i.e. when it merges with Weyl points. It is the p-space analog of a Dirac string terminating on a 
monopole, another analog is given by the Fermi arc in Fig. 1 bottom right . 
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The topological bulk-vortex correspondence exists also for vortices in gapless 
vacua. The topological protection of fermion zero modes is provided by the non- 
trivial topology of three-dimensional Weyl points in the bulk. This bulk-vortex cor- 
respondence II 1001 is illustrated in Fig. |22l In bulk there is a pair of Weyl points 
with opposite topological charges = ±1 in Eq.©. The projections of these Weyl 
points on the direction of the vortex line determine the boundaries of the region 
where the spectrum of fermions bound to the vortex core is exactly zero, E{py ) ~ 0, 
for all within this region. Such flat band was first obtained in Ref. 11241 for the 
noninteracting model. However, due to the topological protection, it is not destroyed 
by interactions. The spectrum of bound states in a singly quantized vortex in ^He-A 
is shown in Fig. |22] The ID flat band terminates at points where the spectrum of 
bound state merges with zeroes in the bulk, i.e. with Weyl points. 



10 Discussion 

The last decades demonstrated that topology becomes a very important tool in 
physics. Topology in momentum space is the main characteristics of the ground 
states of a system at zero temperature (T = 0), in other words it is the characteristics 
of quantum vacua. The gaplessness of fermions in bulk, on the surface or inside the 
vortex core is protected by topology, and thus is not sensitive to the details of the 
microscopic physics (atomic or trans-Planckian). Irrespective of the deformation of 
the parameters of the microscopic theory, the value of the gap (mass) in the energy 
spectrum of these fermions remains strictly zero. This solves the main hierarchy 
problem in particle physics: for fermionic vacua with Fermi points the masses of 
elementary particles are naturally small. 

The vacua, which have nontrivial topology in momentum space, are called the 
topological matter, and the quantum vacuum of Standard Model is the representative 
of the topological matter alongside with topological superfluids and superconduc- 
tors, topological insulators and semi-metals, etc. There is a number of of topological 
invariants in momentum space of different dimensions. They determine universality 
classes of the topological matter and the type of the effective theory which emerges 
at low energy and low temperature. In many cases they also give rise to emergent 
symmetries, including the effective Lorentz invariance and probably all the sym- 
metries of Standard Model, and emergent phenomena such as gauge and gravita- 
tional fields. The symmetry appears to be the secondary factor, which emerges in 
the low-energy corner due to topology, and it is possible that it is topology of the 
quantum vacuum, which is responsible for the properties of the fermionic matter in 
the present low-energy Universe. 

The topological invariants in extended momentum and coordinate space deter- 
mine the bulk-surface and bulk-vortex correspondence. They connect the momen- 
tum space topology in bulk with the real space. These invariants determine the 
fermion zero modes living on the surface of a system or in the core of topologi- 
cal defects (vortices, strings, domain walls, solitons, hedgehogs, etc.). 
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In respect to gravity, the momentum space topology gives some lessons. First, 
in the effective gravity emerging at low energy, the collective variables represent 
the tetrad field and spin connections. In this approach the metric field emerges as 
the composite object of tetrad field, and thus the Einstein-Cartan-Sciama-Kibble 
theory with torsion field is more relevant for the description of gravity (see also 
11125111261 [1271 (1281 ). 

Second, the topology suggests several scenarios of Lorentz invariance violation 
governed by topology. Among them the splitting of Fermi point and development of 
the Dirac points with quadratic and cubic spectrum. The latter leads to the natural 
emergence of the Hof ava-Lifshitz gravity. 
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